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ABSTRACT 

A phenomenological model of X-ray variability of accreting black holes is consid- 
ered, where the variable emission is attributed to multiple active regions/perturbations 
moving radially towards the central black hole. The hard X-rays are produced by in- 
verse Compton upscattering of soft photons coming from reprocessing/thermalization 
of the same hard X-rays. The heating rate of the Comptonizing plasma is assumed 
to scale with the rate of dissipation of gravitational energy while the supply of soft 
photons is assumed to diminish towards the center. Two scenarios are considered: (1) 
an inner hot flow with outer truncated standard accretion disc and (2) an accretion 
disc with an active corona and a thick hot ionized skin. A variant of the model is 
also considered, which is compatible with the currently discussed multi-Lorcntzian 
description of power spectral densities of X-ray lightcurves. 

In the inner hot flow scenario the model can reproduce the observed Fourier fre- 
quency resolved spectra observed in X-ray binaries, in particular the properties of the 
reprocessed component as functions of Fourier frequency. In the accretion disc with 
ionized skin scenario the reduction of soft photons due to the ionized skin is insufficient 
to produce the observed characteristics. 
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1 INTRODUCTION 

X-ray emission from accreting compact objects carries in- 
formation about geometry and physical conditions of the 
accreting plasma. Spectral and timing analyses of the data 
from sources in low/hard state have enabled formulating 
a number of phenomenological scenarios for accretion onto 
black holes. The scenarios belong to two main groups: (1) 
models invoking standard optically thick accretion disc with 
an active corona and (2) models assuming that the optically 
thick disc is truncated at a radius larger than the last stable 
orbit, and is replaced by inner hot, optically thin flow (see 
Done 2002, Czerny 2002 for reviews). Models of both groups 
have some physical foundations. Magnetic dynamos in the 
disc may give rise to the active corona, as suggested by re- 
cent magneto-hydrodynamical simulations (Miller & Stone 
2000), with plasma outflow and/or hot ionized disc skin re- 
sponsible for the observed correlations between spectral pa- 
rameters: the power law slope and amplitude of the reflected 
component (Zdziarski, Lubinski & Smith 1999; Gilfanov, 
Churazov & Revnivtsev 2000); On the other hand, disap- 
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pearance of the cold disc due to evaporation may create the 
inner hot flow (Rozahska & Czerny 2000). Both classes of 
models can explain time average energy spectra, as well as 
correlations between spectral parameters (see discussion of 
models and more references in Zycki 2002; hereafter Z02). 

Specific models of X-ray variability were constructed for 
both geometries (see Poutanen 2001 for review). In the ac- 
tive corona geometry the model of magnetic fiare avalanches 
of Poutanen & Fabian (1999; hereafter PF99) is very suc- 
cessful in reproducing a number of observed properties (e.g. 
power spectra, hard X-ray time lags). However, its simple 
extension to accretion disc with the ionized skin fails to re- 
produce the Fourier-frequency resolved spectra (Revnivtsev, 
Gilfanov & Churazov 1999, 2001) observed from black hole 
binaries (Z02). In the inner hot flow geometry the drifting 
blob model (Bottcher & Liang 1999) reproduces certain ob- 
served characteristics. However, Maccarone, Coppi & Pouta- 
nen (2000) argue that flares in hard (Comptonized) X-rays 
are the primary driver of variability, the soft (thermal) emis- 
sion being merely reprocessed hard X-rays (see also Malzac 
& Jourdain 2000). 

The observed logarithmic energy dependence of the 
hard X-ray time lags was shown to be compatible with a 
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simple propagation model, where X-ray emitting pertur- 
bations move radially towards the central black hole, the 
emitted spectrum hardening as the perturbations approach 
the center (Miyamoto & Kitamoto 1989; Nowak et al. 1999; 
Kotov, Churazov & Gilfanov 2001). Moreover, the observed 
characteristic 'dip' in the energy dependence of time lags 
at the energy of the Fe Ka line can be explained in the 
propagation model, if the line emission from inner disc is 
suppressed (Kotov et al. 2001). 

Therefore in this paper we consider a phenomenologi- 
cal propagation model, where the variable X-ray emission 
is attributed to active regions and/or perturbations moving 
radially towards the central compact object. X-rays are as- 
sumed to be produced by inverse Compton upscattering of 
soft photons coming from a cool, optically thick disc. X-ray 
luminosity is assumed to increase as the emission progresses, 
thus creating a flare of radiation. The Comptonized spec- 
trum is assumed to evolve from softer to harder during the 
flare due to diminishing supply of seed photons. This may be 
due to either the optically thick disc being absent at small 
radii, or the thickness of the ionized disc skin increasing to- 
wards the center. From the simulated event files a number of 
statistics will be computed, which may be compared to avail- 
able data. In particular, the Fourier-frequency resolved spec- 
tral properties are studied: both the continuum slope and Fe 
Kq line strength as functions of Fourier frequency, /, will 
be computed. The observed Fourier-frequency resolved spec- 
tra, calculated from the RXTE data of black hole binaries 
in the low/hard state, show two important characteristics: 
the continuum spectra get harder with increasing /, while 
the amplitude of refiection decreases with / (Revnivtsev et 
al. 1999, 2001). 

We emphasize that the main feature of the considered 
models is the radial propagation of the X-ray emitting struc- 
tures, with the related spectral evolution. This is meant to 
be compared with the models assuming radially localized 
flares (e.g. the model of PF99). However, we do not hy- 
pothesize here about possible correlations between flares (al- 
t hough the idea of flare avalanches will be implemented in 
Sec 3^), hence this model is not meant to explain in detail 
the X-ray power spectra and higher order statistics of the 
variability (Maccarone & Coppi 2002). 

Plan of the paper is as follows: the model is described in 
Sec. ^, results for the hot inner flow geometry are presented 
in Sec. ^ while results for the case of an accretion disc with 
hot ionized skin are presented in Sec. 3.3, In Sec. W a vari- 



ant of the model is considered, which is compatible with 
the multi-Lorentzian description of power spectra (Nowak 
2000). 



2 THE MODEL 

The X-ray emitting structures are assumed to originate 
at a certain radius rout, which will be assumed ~ 50 rg 
(rg = GM/(?\ precise value of rout has little influence on 
our results, since the contribution to emission from large 
distances is insigniflcant). They move towards the center at 
a fraction of the free-fall speed, 



where /? is a model parameter. Since the heating of the 
plasma must ultimately come from the dissipation of the 
gravitational energy, we assume that the plasma heating 
rate scales as the dissipation rate of the latter for a disc- 
like accretion, 



(2) 



where the exponent —2 corresponds to dissipation per ring of 
matter. Here I is the compactness parameter, I 



(^flarc is the characteristic radius of the structure) , and b{r) 
is the boundary term w hich is assumed to have its Newto- 
nian form 6(r) = 1 — 6rg/r (Shakura & Sunyaev 1973). 
Other prescriptions for lh(r) could be used here, since the 
relation between the overall disc accretion flow and the mo- 
tion of structures is obviously uncertain|^ The heating rate 
is normalized to have an assumed maximum value, ih.max- 
Solving the equation of motion (Eq. |l|) gives 



(3) 



where A = | ^^'^ , and the radial positions are expressed 

in units of rg. The duration of a flare is ttrav = ('"out ~ 
Q^/'^yA, since the flnal radial position is assumed 6rg. The 
resulting flare, /h(i), is characterized by a slow rise and a 
sudden decrease (Fig. |l^) in agreement with the constraints 
derived by Maccarone et al. (2000). 

The soft photons are assumed to come from reprocess- 
ing and thermalization of the hard X-rays, so that the feed- 
back loop is realized as needed to explain the R-V correla- 
tion (Zdziarski et al. 1999; Gilfanov et al. 2000). Two geo- 
metrical scenarios may be considered here: 

(i) an inner hot flow partially overlapping a cold, optically 
thick disc disrupted at certain radius, rtr, (see e.g. Pouta- 
nen, Krolik & Ryde, 1997, for arguments for the overlap). 
The luminosity of soft photons crossing an active region is 
therefore parametrized as 



for r > rtr 
for r < rtr. 



(4) 



V = Pvf[ = P 



2GM 



(1) 



where A'^b ~ 0.5, 

(ii) a standard optically thick accretion disc with an ac- 
tive corona and a hot ionized skin on top of the disc (Nayak- 
shin & Dove 2001). The flux of soft photons crossing an 
active region is computed using the method of Zycki & 
Rozahska (2001) and Z02, for an assumed height of the ac- 
tive region above the disc, as it travels towards the center. 

It should be noted here that for a truncated disc, wholly 
replaced by a hot flow, and for a compact emission region, 
the resulting dependence Zs(r)/Zh(r) is much steeper than 
the assumed r^. That is, once the active region is within the 

^ One subtle geometrical effect here is the radial dependence of 
the size of the emiter, -Rfjare- Actually, it is the luminosity from 
each emitter that can be expected to scale as r~^. Then, eq. ^ is 
equivalent to an assumption -Rflare = const. If -Rflare varies with r, 
radial dependencies of L and will be different, but the observed 
quantity is L and this should scale as . Precise values of /h(r) 
is not important in our paper because of approximations used 
to compute the spectra (eq. In particular, pair production 
process es a, re neglected. This problem will be more apparent in 
Section 3.3 
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Figure 1. Time variations of parameters during a flare in the 
truncated disc scenario, (a) luminosity of a flare, (b) spectral 
index, F, and amplitude of reflection, R. Parameters: rout = 60 Tg, 
13 = 0.02 so that the flare duration is ttrav = 0.5 sec, rtr = 30rg. 



central hole, it can intercept only a tiny fraction of the ther- 
malized emission, if its size is comparable or smaller than the 
transition radius. This would produce much harder spectra 
than observed, F < 1. Therefore we choose here a weaker 
dependence, although the particular value of the exponent 
is chosen rather arbitrarily and no attempt is made here to 
tune it to the observed data. The slow dependence may be 
interpreted as the optically thick plasma slowly disappearing 
(disc disrupted into small clouds; covering factor gradually 
decreasing from 1 to 0) with decreasing r, or the emission 
coming from perturbations propagating in the entire volume 
of the hot plasma. 

The emitting structures are assumed to originate at 
Tout > rtr at a mean rate of A per second. Time intervals 
between the launch of successive flares are generated from 
the Xexp{ — Xt) distribution, as appropriate for a Poissonian 
process. Each flare is followed until it reaches the marginally 
stable orbit at rms = 6rg. At each time step the primary 
Comptonized component and its reprocessed component are 
computed. The truncation radius is assumed rtr ~ 30 rg, so 
that the time average spectral slope of the Comptonized 
continuum is F = 1.6-1.7. 

The Comptonized spectrum is computed using the code 
thComp (Zdziarski, Johnson & Magdziarz 1996), solving 
the Kompaneets equation. Computations are parametrized 
by the photon spectral index, which we compute from 
F = 2.33(Zh/Zs)"^''® (Beloborodov 1999a,b), and electron 
temperature kTe/ {rne<?) computed using formulae from Be- 
loborodov (1999b). Plasma optical depth is assumed Tcs = 
1.8 (PF99). 

The reprocessing medium is assumed to be 'cold', (only 
hydrogen and helium ionized) with elements abundances of 
Morrison & McCammon (1983). The spectrum of the Comp- 
ton reflected continuum is computed from the simple for- 
mula of Lightman & White (1988), 



^abs 



^abs "i" ^c. 



(5) 



where Kabs(i5) and Kcs{E) are the photo- absorption and elec- 
tron scattering opacities, respectively. The above formula is 
multiplied by a simple exponential cutoff to mimic the high 
energy cutoff. This is sufficient for our purposes since only 
the spectra up to 30 keV will be considered. The Fe Ka 
line is added to the reflected continuum, with equivalent 



width (EW) as a function of F, following computations of 
Zycki & Czerny (1994; see also Z02). Importantly, the in- 
stantaneous amplitude of the reprocessed component is as- 
sumed to depend on radius of emission. 



R{t) 



for r > rtr 
for r < rtr, 



(6) 



i.e. it linearly decreases from 1 at r = rtr to at r = rms. 
Again, the particular form of the function R{r) for r < rtr 
is chosen rather arbitrarily, the general condition being that 
R{r) decreases inwards. Time evolution of F and R during 
a flare is plotted in Fig |l|b. 

Thus created sequence of spectra is subject to standard 
analysis in the time and Fourier domains (see e.g. van der 
Klis 1995; Nowak et al. 1999; Poutanen 2001). In partic- 
ular, the Fourier frequency resolved (/-resolved hereafter) 
spectra are computed according to the prescription given by 
Revnivtsev et al. (1999; see also Z02). The normalized power 
spectral density (PSD) is computed as 



J =0,...,iV-l 



(7) 

(8) 



where the discrete frequencies fj = j/T for j — 
-N/2, . .., N/2; T is the total light curve time, C is the 
mean count rate, while Ck are the number of counts in fc-th 
time bin. The /-resolved spectrum at energy Ei and Fourier 
frequency fj is then deflned as 



S{E,J,)=a^Pi{fj)Af, 



-A/, 



(9) 



All computations are done assuming the central black 
hole mass M = 10 . 



3 RESULTS 

3.1 The case of a single value of velocity of the 
active regions 

In the simple case of a single value of /3 we do not expect 
the observed PSD to be reproduced, since the behaviour 
generally requires a distribution of parameters and/or cor- 
relations between flares (Lehto 1989, Lochner, Swank & 
Szymkowiak 1991). It is nevertheless instructive to study 
other observables from the model. The velocity proportion- 
ality coefficient was adopted P = 0.02 and rout ~ 60 rg, so 
that the duration of a flare is ftrav ~ 0.5 s, and the PSD 
peaks at 0.5-1 Hz. Fig. |^ shows the /-resolved spectra and 
time lags as functions of energy. The model predictions do 
not match well the observed trends. The time lags as a func- 
tion of energy are independent of Fourier frequency, at least 
for / < 1 Hz, that is below the peak frequency of PSD. This 
contrasts the observed 5t{E), which shows longer lags at a 
given E, for lower / (Nowak et al. 1999). The logarithmic de- 
pendence of energy is true for any flare proflle, if the change 
of spectral index during the flare, AF, is small compared to 
time average F (Kotov et al. 2001). In the considered case 
AF ~ 0.8, so it is not small compared to (F) « 1.7, but 
the logarithmic dependence St{E) oc logiJ is still obeyed. 
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Figure 2. Results of simulatio ns f or a single value of active regions velocity. Parameters: /3 = 0.02, rout = 60 rg, so that the duration 
of a flare is ttrav ~ 0.5 s (Sec. 3.1). (a) The /-resolved snectra, divided by a power law, Ne oc E~^'^ , and normalized to 1 at 2 keV 
(labels are middle frequencies in A/j in Hz as in equation Bj) The spectra harden initially with / up to / ~ 10 Hz, then soften, with the 
corresponding increase of the reflection amplitude, (b) Equivalent width of the Fe Ka line showing dependence on / matching that of 
the /-spectra, (c) Time lags as functions of energy, for three values of Fourier frequency (Kotov et al. 2001; also fig. 12 in Nowak et al. 
1999). 
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Figure 3. Results of simulations for avalanches of flares. Curves are described in the same way as in Fig. y. Duration of the flares covers 
the range 0.03 s to 1.5 s (/3 = 0.35-0.007, respectively). The /-resolved spectra and EW of Ka line are now monotonic functions of 
Fourier frequency /, similarly to the observed dependencies (Revnivtsev et al. 1999). The 5t{E) dependence, panel (c), reproduces the 
observed dependencies well (compare with fig. 12 in Nowak et al. 1999). 
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Figure 4. Results of simulations for a simple distribution of velocities of emission regions (no correlations between fiares, i.e. no 
avalanches). Curves are described in the same way as in Fig. ^ The /-resolved spectra and EW of Ka line are very similar to the 
avalanche case (Fig. but the time lags St{E) are much longer and longer than those observed in Cyg X-1. Longer individual flares 
are necessary to generate variability power at frequencies down to 0.1 Hz. 
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at least for higher / > 1 Hz, and energies not exceeding 
~ 30 keV. This is consistent with analytical solution, which 
can be obtained for a simple exponential profile and a linear 
dependence of T{t). Assuming 

S{E,t) =exp{t/t,) E-'^f+''\ t<0, (10) 

where F/ is the spectral index at the peak of a flare, at t = 0, 
the Fourier transform can be computed. 

The phase lags are 

<^(/) = arg[S(i5i,/)S*(i52,/)], (12) 

where 

tan<^(/) = 

-a«gln(Ei/-E2)2'r/ , ^ 

1 + atr ln(EiB2) + i2tr InBi In B2 +<r (2'r/)2 ' ^ ' 

From this it can be seen that for small (j>, the time lags St{E) 
are independent of /, 5t = (j)/{2nf) « tan (j) / {2n f) ~ const. 

The dip in St{E) at the energy of the Fe Kq line, seen 
in the real data, is reproduced by the model, confirming the 
result of Kotov et al. (2001). 

The observed spectral slope of the comptonized contin- 
uum and the equivalent width of the Fe Kq line are mono- 
tonic functions of /, at least above ~ 0.5 Hz (Revnivtsev et 
al. 1999). In contrast to these, the model /-resolved spec- 
tra show complex dependence on /. In particular the EW 
of the Ka line first decreases above ~ 1 Hz, then increases 
again above ~ 10 Hz. This means that in the Fourier de- 
composition of the flare profile the highest frequencies form 
the first part of the rising profile, which, in real space, is 
formed above the accretion disc (i.e. rti < r < rout). The 
spectra from that part are soft, and the amplitude of re- 
fiection is close to 1. However, we note that the /-resolved 
spectra are dependent on the radial profile of the heating 
function (Eq. Q). Experimenting with the lh{r) function we 
find that a more peaked lh{r) (e.g. oc r 7 > 3) produces 
the /-resolved spectra monotonically hardening with / and 
EW(/) monotonically decreasing. 

3.2 Distribution of parameters 

The full model obviously has to explain all observable char- 
acteristics of the X-ray emission. The complex but universal 
PSD of X-ray emission in low/hard state points out towards 
a universal mechanism of generating the complexity. One 
such mechanism was suggested by Stern & Svensson (1996) 
in the context of gamma-ray bursts and implemented by 
PF99 for accreting black holes. Their idea is that each flare 
can stimulate one or more flares to produce an 'avalanche' 
of many flares. Parameters of avalanches can be adjusted so 
that the power spectrum can have the observed slope of — 1 
(i.e. P{f) oc f~^) in a given frequency band. 

The idea of avalanches of flares may be easily imple- 
mented in the present model. Difl'erent time-scales of indi- 
vidual flares correspond to difl'erent coefficients /3, and the 
stimulated flares are assumed to originate from the same 
radius as the stimulating ones. Specifically, we assume that 
(3 follows a power law distribution P{(3) oc between f3i 
and (32 ■ There are A spontaneous fiares per second, and each 



flare has a probability /i of stimulating a secondary flare, the 
latter being delayed by tdeiay = aitrav/7, where ttrav is the 
duration of the stimulating flare, while a is drawn from a 
Poisson distribution of mean qq (the factor 7 approximately 
translates ftrav to flare time scale used by PF99; see there for 
detailed description of the model) . The adopted values of the 
parameters are such that they give good description of PSD 
from black hole binaries: k = 0.5, /3i = 0.007, /32 = 0.35, 
A = 90, fi — 0.7, ao = 5. The perturbations start from 
fout ~ 60 rg, so the duration of the flares cover the range 
from 0.03 to 1.5 sec. 

Results for this model are presented in Figure ^ They 
match the observed dependencies reasonably well. Time lags 
as functions of energy are in good quantitative agreement 
with the observed values in the whole range of computed 
/ (compare with Cyg X-1 data in Nowak et al. 1999). In- 
terestingly, by superposing flares of difl'erent durations, the 
rising trend in EW(/) for / > 10 Hz (see Fig. |^ was com- 
pensated for, and the equivalent width is now monotonically 
decreasing with /. Correspondingly, the continuum hardens 
monotonically with /. Both trends are as observed in black 
hole binaries (Revnivtsev et al. 1999, 2001). This means that 
not only are higher frequencies formed from shorter flares, 
but within the flares the dominant contribution to a given 
/ moves towards the peak of the flare, as / increases. 

For comparison. Fig. ^ shows results obtained assum- 
ing a simple distribution of velocities of the emission struc- 
tures (no correlations between flares). Coefficients (3 was 
drawn from a uniform distribution between /3i = 0.003 
and 132 = 0.1 (fiare durations 0.1 s to 3.5 s). Importantly, 
the lower value /3i is now much smaller than that in the 
avalanche prescription, to give fiat P{f) down to « 0.1 Hz. 
As a result, the time lags are now longer than those for 
avalanches, clearly longer than those observed. 

3.3 Accretion disc with hot ionized skin 

The presented propagation model can also be formulated in 
the geometry of an accretion disc extending down to the 
last stable orbit, with the X-ray emitting structures mov- 
ing above the disc towards the center. As a result of ir- 
radiation by hard X-rays, the hot ionized skin will form 
on the disc (Rozanska & Czerny 1996; Nayakshin & Dove 
2001). Making the same assumption about increasing lu- 
minosity of the moving flare as in Section ^ (equation 
the result is obtained that the thickness of the hot skin, 
Thot, increases towards the center (Nayakshin 2000; Zycki & 
Rozanska 2001). This reduces the effectiveness of reprocess- 
ing/thermalization as the flare progresses. This in turn leads 
to diminishing supply of soft photons for Comptonization 
and hardening of the continuum as a consequence. Simul- 
taneously, the relative amplitude of the "cold" reprocessed 
component decreases, since fewer photons penetrate to the 
cold disc beneath the hot skin (e.g. Done & Nayakshin 2001; 
Zycki & Rozanska 2001). Qualitatively, the situation is thus 
analogous to the disrupted disc scenario, and could repro- 
duce observed trends of F and R with Fourier frequency. 
In this Section a quantitative analysis of the model is per- 
formed. 

Basic model element is the spectral evolution during 
a single flare. Here the radial dependence of the disc il- 
luminating flux, -Firr('"), is crucial since it determines the 
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Figure 5. Time var iati ons of parameters during a flare in the hot 
skin scenario (Sec. [3.3| ). Left panel: luminosity of a flare (solid 
curve) and the thickness of the hot skin (dashed curve), right 
panel: spectral index, F, and amplitude of reflection, R. Note the 
relatively small range of F and R, due to small variation of rjjot. 
The spectral evolution during a flare is insufficient to produce 
correct hard X-ray time lags or /-resolved spectra. 



Thot('') dependence. The illuminating flux can be written as 
Fin ~ {t]/tv)L/ {R^^j.^+SH'^), where rj « 0.5 is the anisotropy 
factor, -Rflaro is the radius of a flare and H is its height above 
the disc (see PF99). Since, as argued above (equation 



L{r) oc 



we would obtain -Firr('') oc r if i?flarc = const 



and H = const. However, for Thot to have a significant radial 
dependence, _Firr must be a stronger function of r. For exam- 
ple, for illumination due to local dissipation, Firri'r) oc r~^. 
To achieve this, we assume here Raa.icir) oc H{r) oc r^^^. 
That is, the flare shrinks and its height decreases as it ap- 
proaches the center. The condition -Rfiarc oc // is required in 
order to reproduce the T-R correlation, since it ensures that 
the geometry of the flare is constant and the only varying 
factor is the thickness of the hot skin. The proportionality 
coefficient in the -Rflaro = Ror^^^ relation can be com- 
puted if the maximum heating compactness, total source 
luminosity and fiare occurence rate are known (see eq. 7 in 
Z02). 

Energy spectra during a flare are computed as in Z02. 
The thickness of the hot skin, Thot, is computed as a func- 
tion of Fin and r using the fitting formula (eq. 9) from that 
paper, while the amount of soft photons from thermalization 
re-entering an active region and the effective amplitude of 
reflection were computed from eqs. 3, 11 and 10 in Z02, re- 
spectively, following PF99 and Poutanen (2002). The feed- 
back parameter Do = 0.5 (eq. 3 in Z02 and PF99)|. The 
height H is assumed to be 0, so that in the absence of the 
hot skin F ~ 2. It should be emphasized that in the present 
scenario the time/radial dependencies of Ih/L and R can be 
predicted from physical considerations. Thus, eq. ^and eq. ^ 
are not used here. 

An example of a flare computed in this model is shown 
in Fig H The spectrum does evolve from soft to hard, but 

We note that the parametrization adopted by PF99 hides the 
fact that the required Dg is actually incompatible with an emis- 
sion region whose horizontal and vertical extent are comparable, 
if the emission is isotropic. According to PF99 for Dg = 0.5 one 
half of the emitted X— rays return to the emission region if H = 0. 
This in reality implies the horizontal extent much larger than 
vertical one (slab-like geometry). 



the change of F is very small, clearly insufficient to produce 
any difference in /-resolved spectra. A sequence of spectra 
was simulated, assuming a distribution of velocities of active 
regions and other parameters as in pulse avalanche prescrip- 
tion (Sec. 3.2). Computing the /-resolved spectra we confirm 
that they all have the same slope, F « 2, while the EW of 
the Ka line is approximately constant at ~ 55 eV. Thus, the 
hot skin disc model in the simple version considered here is 
not able to reproduce the /-resolved spectra. 

The irradiation flux would have to be a very strong 
function of r for the hot skin thickness to increase enough 
to produce the required hard spectrum (F ~ 1.6-1.7 at the 
peak of the flare). This would mean a signiflcant decrease of 
the size of the emitting region as it approaches the center, if 
the luminosity L{r) oc dependence is kept. That in turn 
would increase the compactness parameter, Zh, as well as the 
temperature of the soft thermalized radiation. The former 
cannot be uniquely constrained by the data, if the plasma 
is thermal (Maciolek-Niedzwiecki, Zdziarski & Coppi 1995), 
but the latter is observed to be no higher than To ~ 0.3 keV 
(e.g. Di Salvo et al. 2001). It remains to be seen whether a 
physically plausible model can be constructed fulfilling all 
the constraints. 



4 A MODEL COMPATIBLE WITH 

MULTI-LORENTZIAN DECOMPOSITION 
OF POWER SPECTRA 

X-ray power spectra from accreting compact objects can be 
interpreted in terms of a (small) number of relatively nar- 
row Fourier components (Nowak 2000; Belloni, Psaltis & van 
der Klis 2002). Higher quality data from RXTE clearly re- 
veal that PSD, which were previously considered a roughly 
featureless power law with two or three breaks, consist of 
a number of peaks. This interpretation is strengthened by 
the correlations between frequencies of the components (Bel- 
loni et al. 2002), reducing the complexity of broad band 
variability to possibly just one parameter, at least within 
a given spectral state. The parameter seems to be related 
to the geometry of accretion fiow, as evidenced by correla- 
tions between the characteristic frequencies and energy spec- 
tral properties (spectral slope, amplitude of the reprocessed 
component; Revnivtsev et al. 2001; Gilfanov et al. 2000). 
The obvious interpretation of the geometrical factor is the 
variable transition radius between the outer optically thick 
standard accretion disc and inner optically thin hot flow 
(Revnivtsev et al. 2001; Psahis & Norman 2002). 

The narrow components of PSD are usually modeled as 
Lorentz profiles. If the Lorentz profiles are taken seriously, 
the emission may be envisioned to consists of a collection of 
exponentially damped (or forced i.e. rising) harmonic oscil- 
lators. Obviously, the oscillations may be reduced to simple 
exponential shots, if the quality factor of the Lorentzians, 
Q ~ A///o (/o is the resonant frequency while A/ is the 
half- width of the Lorentzian), is low enough, Q < 1. In the 
interpretation of Psaltis & Norman (2002), the frequencies 
would be connected to characteristic frequencies at the disc 
truncation radius: orbital, periastron-precession and nodal- 
precession (Stella & Vietri 1998). However, it is the hard 
X-ray emission that is actually variable, and the most likely 
reason for the variability is a variable energy dissipation 
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in the hot plasma, rather than variable input of soft pho- 
tons (Maccarone et al. 2000; Malzac & Jourdain 2000), even 
though it is far from obvious how to make a connection be- 
tween the disc frequencies and the energy dissipation in the 
hot plcisma. Furthermore, hard X-ray time lags can be in- 
terpreted as a result of spectral evolution during individual 
emission events (shots/flares in the usual shot noise models; 
PF99; Kotov et al. 2001). 

In this Section we consider a descriptive model com- 
patible with the above ideas. Adopting the form of the 
Lorentzian function as in Nowak (2000), 



N^Qfo 



fS+QHf~for 



(14) 



where A*', Q and /o are amplitude, quality factor and reso- 
nant frequency of the oscillations, we simulate the shots of 
the following temporal profile: 



A{t) = Ao exp (^-^^ cos(27r/ot -I- 0o), t <0. 



(15) 



Here, the normalization Ao — NQ^^^, the rise time tr ~ 
2Q/(27r/o) and 0o is a randomly chosen initial phase. Pa- 
rameters of the Lorentzians are taken from fits to Cyg X-1 
0-4 keV light curves (Table 2 in Nowak 2000) . We ignore the 
lowest amplitude component QPO2, which leaves us with 5 
components. 

The crucial assumption we make here is that the shots 
correspond to perturbations traveling through the plasma, 
from the truncation radius, rti, towards the center, at a 
speed which may, but does not need to be, the same for all 
components. In actual computations the speed is assumed 
constant in time, but assuming v oc vs (as in previous Sec- 
tions) would not change our conclusions. The heating rate 
of the plasma is assumed to be proportional to the energy 
of the oscillations, i.e. lh{t) oc [j4(t)]^, as in usual oscillatory 
motion (the square of a damped/forced oscillator is again a 
damped/forced oscillator with a Lorentzian PSD). We em- 
phasize that the heating is now not described by eq. ^, since 
it is determined by A{t), but formulae for ls{t) and R{t) 
(eq. 1^ and eq. ^ are adopted. Each of the 5 components is 
chosen with the same probability, while the overall frequency 
of shots is adjusted to give the r.m.s. variability « 30 per- 
cent. Computations of spectra and subsequent Fourier anal- 
ysis proceeds as in previous models. 

In the particular case considered here the longest ti ~ 
1.2 sec for the lowest /o component, and ti is progres- 
sively shorter for the higher- /o components. The propaga- 
tion speed, V, is chosen such that the travel time, ttrav ~ 
{rtT — rf)/v is at least 5ir. This means that a shot is launched 
from rtr at f = — Sir, where is the corresponding rise time- 
scale (the beginning of an exponential shot is obviously not 
well defined). The final minimum radius, rf , reached at t = 0, 
is determined by v. 

It is easy to predict qualitatively the results for /- 
resolved spectra, if the propagation speed is the same for 
all components. The higher the /o (the shorter the U), 
the shorter the distance traveled. Shorter distance means 
softer energy spectrum, since spectra produced closer to the 
truncation radius are, by construction, softer. Thus the /- 
resolved spectra can be expected to become softer with in- 
creasing /, opposite to those observed. Thus, obviously, a 
necessary condition for this model to reproduce the observed 



/-spectra is to assume that the speed of the perturbations 
increases with fo- Then, the higher the frequency /, the far- 
ther from the cold disc the emission takes place, producing 
harder spectra. Adjusting the speeds accordingly we indeed 
obtain the required dependence of /-spectra (including EW 
of the Ka line) on Fourier frequency. 



5 DISCUSSION 

We have considered a propagation model of X-ray variability 
in application to high energy emission from accreting com- 
pact objects. The motivation for our considerations comes 
from two directions: (1) problems faced by the models of 
localized flares (PF99) in explaining the /-resolved spectra 
(Z02), and (2) demonstration by Kotov et al. (2001) that 
a propagation model with spectral evolution can naturally 
explain the hard X-ray time lags. The model can be formu- 
lated in both considered geometries: truncated standard disc 
with inner hot flow and a standard disc with an active corona 
and the hot ionized skin. It should be emphasized, however, 
that the main parameters entering the spectral computa- 
tions, lh/ls{t) and R{t) can be predicted only in the latter 
model, while they have to be assumed rather arbitrarily in 
the former. It is the greater predictive power of the hot skin 
scenario what enables to find unique model predictions and 
demonstrate the problems which the model faces. 

Similar computations of variability could be performed 
in the context of the plasma outfiow model of Beloborodov 
(1999a,b). Here the main parameter is the outfiow speed. 
Qualitatively, the requirement for the model to work is that 
the speed increases inwards. A good physical model predict- 
ing the speed as a function of distance. X-ray fiux and other 
parameters is necessary, in order to make quantitative pre- 
dictions. 

Let us concentrate now on the geometry of an inner 
hot X-ray producing fiow with an outer standard optically 
thick disc. This is one of the possibilities of the geometri- 
cal arrangement of accreting multi-phase plasma. It natu- 
rally accounts for hard spectra of black hole X-ray bina- 
ries, small amplitude of the reprocessed component, and ob- 
served correlations between these spectral parameters (Done 
2002 and references therein). Emission in the form of flares 
(shots) associated with active regions and/or perturbations 
traveling in the hot plasma might then be responsible for 
observed time and spectral variability, considering that X- 
ray variability is generally of stochastic character (Czerny 

6 Lehto 1997). Such an extension of models of steady emis- 
sion from hot flows (e.g. Esin, McClintock & Narayan 1997) 
is necessary, considering the many observable characteris- 
tics related to variability. The situation is fully analogous to 
models invoking active regions above an accretion disc (e.g. 
Stern et al. 1995). These models could explain the steady 
Comptonized emission by basically adjusting one parame- 
ter, the ratio of heating to cooling compactnesses. When 
supplemented with the idea of plasma outfiow (Beloborodov 
1999a,b) or a hot ionized skin (Nayakshin & Dove 2001), 
they cloud explain the diversity of spectra and correlations 
between spectral parameters. However, in order to explain 
the complex time variability, a number of additional features 
had to be added: correlations between fiares (avalanches) to 
explain the PSD and spectral evolution during a fiare as a 
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result of, e.g., rising the active region above the disc (PF99). 
Similar developments can be expected in the hot flow sce- 
nario. 

Precise geometry of the X-ray emitting structures in 
the hot flow is obviously uncertain, but certain constraints 
can be obtained. For example, compact active regions (size 
much smaller than the inner disc radius) moving through hot 
plasma totally devoid of optically thick plasma (i.e. the cold 
disc sharply disrupted at certain radius), would intercept so 
few soft photons at r < rtr that the continuum slope would 
be much harder than observed. One possibility would then 
be that the active regions are compact but the disappear- 
ance of the optically thick phase is gradual, moving from rtr 
inwards. This could mean that small optically thick clouds 
form inside the hot plasma at r < rtr and their fllling factor 
gradually decreases inwards from ~ 1 to 0. Detailed sim- 
ulations would be necessary to predict spectra from such a 
configuration (e.g. Malzac & Celotti 2002). At a rather spec- 
ulative level one might point out two possibilities of forma- 
tion of such compact traveling active regions. The accreting 
plasma might form multiple shocks where the dissipation of 
energy would preferentially take place, or the active regions 
might have something to do with the dissipation of magnetic 
energy inside the hot flow, which dissipation is necessary to 
maintain the equipartition between the magnetic and ther- 
mal energy (Bisnovatyi-Kogan & Lovelace 2000). 

Another possibility for the geometry is to assume that 
the disc is rather sharply truncated but the X-ray emis- 
sion comes from global coherent perturbations propagating 
through the entire volume of the hot plasma. Then the de- 
crease of Is/W would roughly correspond to the decrease of 
the solid angle subtended by the cold disc, as the perturba- 
tion approaches the center. 

Temporal profiles of the flares considered in this paper 
were rather simple, therefore to produce a broad band PSD 
a distribution of parameters (e.g. flare timescale) was re- 
quired. One might envision more complex profiles based on, 
e.g., recent considerations of Maccarone & Coppi (2002). 
They suggest correlations between launch times of short 
flares (tr ~ 0.5 sec), such that the short flares flU in a spe- 
ciflc temporal envelope function. We note that their flares 
do not explain the whole range of time-scales, since the PSD 
above ~ 1 Hz is under-predicted. Even shorter flares would 
be necessary to generate enough power at sub-second time- 
scales, perhaps fllling in the short flares considered by Mac- 
carone & Coppi (2002) as envelopes. The result might then 
be equivalent to a hierarchical scenario proposed by Uttley 
& McHardy (2001), where a long time-scales perturbation 
propagate from outside, breaking down into smaller, faster 
structures in a fractal-like manner. 



ACKNOWLEDGMENTS 

This work was partly supported by grant no. 2P03D01718 of 
the Polish State Committee for Scientific Research (KBN). 



REFERENCES 

Belloni T., Psaltis D., van der Klis M., 2002, ApJ, 572, 392 
Beloborodov A. M., 1999a, ApJ, 510, L123 



Beloborodov A. M., 1999b, in Poutanen J., Svensson R. eds, ASP 

Conf. Ser. 161. High Energy Processes in Accreting Black 

Holes, 295, ( astro-ph/990110^ ) 
Bisnovatyi-Kogan G. S., Lovelace R. V. E., 2000, ApJ, 529, 978 
Bottcher M., Liang E. P., 1999, ApJ, 511, L37 
Czerny B., 2002, in Active Galactic Nuclei: from Central Engine 

to Host Galaxy, ASP Conference Series, Eds. S. Collin, F. 

Combes and L Shlosman, in press 
Czerny B., Lehto H. J., 1997, MNRAS, 285, 365 
Di Salvo T., Done C, Zycki P. T., Burderi L., Robba N. R., 2001, 

ApJ, 547, 1024 

Done C, 2002 , Philosophical Tr ansactions of the Royal Society, 

360, 1967 ( |astro-ph/0203246| ) 
Done C, Nayakshin S., 2001, MNRAS, 328, 616 
Esin A. A., McClintock J. E., Narayan R., 1997, ApJ, 489, 865 
Gilfanov M., Churazov E., Revnivtsev M., 2000, Proc. 5th Sino- 

German Workshop on Astrophysics, SGSC Conference Ser., 

Vol. 1. China Scie nce & Technology Press, Beijing, p. 114 

( astro-ph/000241E| ) 
Kotov O., Churazov E., Gilfanov M., 2001, MNRAS, 327, 799 
Lehto H. J., 1989, in The 23rd ESLAB Symposium on Two Topics 

in X-Ray Astronomy. Volume 1: X Ray Binaries, ESA, p. 499 
Lightman A. P., White T. R., 1988, ApJ, 335, L57 
Lochner J. C, Swank J. H., Szymkowiak A. E., 1991, ApJ, 376, 

295 

Maccarone T 



ph/0204153) 



Coppi P. S., 2002, MNRAS, in press (astro 



Maccarone T. J., Coppi P. S., Poutanen J., 2000, ApJ, 537, L107 
Maciolek-Niedzwiecki A., Zdziarski A. A., Coppi P. S., 1995, MN- 
RAS, 276, 273 
Malzac J., Celotti A., 2002, MNRAS, 335, 23 
Malzac J., Jourdain E., 2000, A&A, 359, 843 
Miller K. A., Stone J. M, 2000, ApJ, 534, 398 
Kitamoto S., 1989, Nat, 342, 773 
McCammon D., 1983, ApJ, 270, 119 
2000, ApJ, 534, 718 
Dove J. B., 2001, ApJ, 560, 885 
Nowak M. A., 2000, MNRAS, 318, 361 

Nowak M. A., Wilms J., Vaughan B. A., Dove J. B, Begelman M. 

C, 1999, ApJ, 515, 726 

Poutanen J., 2001, AdSpR, 28, 267 ( astro-ph/0102325 ) 
Poutanen J., 2002, MNRAS, 332, 257 

Poutanen J., Fabian A. C, 1999, MNRAS, 306, L31 (PF99) 
Poutanen J., Krolik J. H., Ryde F., 1997, M NRAS. 292. L21 



Miyamoto S. 
Morrison R., 
Nayakshin S., 
Nayakshin S., 



Psaltis D., Norman C, 2002, ApJ, in press ( |astro-ph/0001391 ) 
Revnivtsev M., Gilfanov M., Churazov E., 1999, A&A, 347, L23 
Revnivtsev M., Gilfanov M., Churazov E., 2001, A&A, 380, 520 
Rozanska A., Czerny B., 1996, Acta Astron., 46, 233 
Rozariska A., Czerny B., 2000, A&A, 360, 1170 
Shakura N. I., Sunyaev R. A., 1973, A&A, 24, 337 
Stella L., Vietri M., 1998, ApJ, 492, L59 
Stern B. E., Svensson R., 1996, ApJ, 469, L109 
Stern B. E., Poutanen J., Svensson R., Sikora M., Begelman M. 

C, 1995, ApJ, 449, L13 
Uttley P., McHardy I. M., 2001, MNRAS, 323, L26 
van der Klis M., 1995, in Lewin W. H. G., van Paradijs J., van 
den Heuvel E. P. J., eds. X-ray binaries, Cambridge Univ. 
Pres, Cambridge, p. 252 
Zdziarski A. A., Johnson W. N., Magdziarz P., 1996 ApJ, 283, 
193 

Zdziarski A. A., Lubinski P., Smith D. A., 1999, MNRAS, 303, 
Lll 

Zycki P. T., 2002, MNRAS, 333, 800 (Z02) 
Zycki P. T., Czerny B., 1994, MNRAS, 266, 653 
Zycki P. T., Rozanska A., 2001, MNRAS, 325, 197 



© 0000 RAS, MNRAS 000, 000-000 



